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(-H We report on recent developments in the public computer code HiggsBounds, which confronts 

arbitrary Higgs sector predictions with 95% C.L. exclusion limits from Higgs searches at the LEP, 
Tevatron and LHC experiments. We discuss in detail the performance of the Standard Model 
(SM) likeness test as implemented in the latest version HiggsBounds-3 .8.0, whose outcome 
decides whether a search for a SM Higgs boson can be applied to a model beyond the SM. Fur- 
thermore, we give a preview of features in the upcoming version HiggsBounds-4 .0.0 and 
■ the new program HiggsSignals, which performs a test of Higgs sector predictions against 

the signal rate and mass measurements from Higgs boson analyses at the Tevatron and LHC. This 
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is illustrated with an example where the heavier CP-even Higgs boson of the Minimal Super- 
symmetric Standard Model (MSSM) is considered as an explanation of the LHC Higgs signal at 
~ 126 GeV. 
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1. Introduction to HiggsBounds 

The search for the Higgs boson at coUider experiments is essential for unravelhng the mech- 
anism of electroweak symmetry breaking. These searches have been carried out over the last two 
decades at the LEP and Tevatron experiments, and for the last two years also at the LHC experi- 
ments. Negative search results^ are used to derive 95% C.L. limits on the Higgs signal rate. These 
results are provided either as fairly model-independent limits on the cross section of a single sig- 
nal topology {e.g. the LEP process e^e^ — )• /i,Z — )• bbZ or the charged Higgs process in top quark 
production, t — )■ H^b — )• TV^b) or as a limit on combined cross sections for a specific model, in 
particular the Standard Model (SM). In the latter case, the discovery potential is maximized by 
the combination of several production and/or decay modes, at the cost of a more model-dependent 
cross section limit. 

The limit on the Higgs boson signal rate is derived by rejecting the Higgs signal plus back- 
ground hypothesis. An example of a model-dependent exclusion limit is given in Fig. 1 for the 
ATLAS search for a SM Higgs boson decaying to two photons [1]. The figure presents both the 
observed limit and the expected limit on a universal scale factor \i = o j asM of the combined SM 
Higgs boson signal rate as a function of the assumed Higgs mass mu. The expected limit is derived 
from Monte Carlo simulation where the data is assumed to be identical to the background expecta- 
tion. In this example, the considered Higgs production modes comprise of the gluon-gluon fusion 
(ggf) and vector boson fusion (VBF) process as well as the associated Higgs production with a 
vector boson (HW , HZ) and with a top quark pair {Htt). Without the knowledge of the individual 
efficiencies of these signal topologies, the limit applies only to models where the five production 
modes ai^e approximately in the same relations as in the SM. This is discussed in detail in Section 2. 

The public computer code HiggsBounds [2, 3, 4] confronts arbitrary Higgs sector predic- 
tions with the most recent 95% C.L. exclusion limits from neutral and charged Higgs searches from 
the LEP, Tevatron and LHC experiments. In order to preserve the 95% C.L. limit interpretation. 
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Figure 1: 95% C.L. exclusion limit from the ATLAS Standard Model H ^ yy search [1]. The limit is given 
on the signal strength modifier, \i = ct/ctsm, which universally scales the Standard Model Higgs boson 
signal rates of the considered production and decay modes. 



'We discuss the recent discovery in the Higgs searches at the LHC in Sect. 3. 
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we follow a well-defined statistical procedure: First, the most sensitive analysis {i.e. the analysis 
with the best expected exclusion limit for the particular parameter point of the investigated model) 
is determined out of all available analyses. Then, the parameter point is tested against the observed 
limit of this particular analysis only. 

As input the code requires the number of neutral and charged Higgs bosons hi, their masses 
m/,., total widths Txox{hi), decay branching ratios BR(/j; —;•...), (normalized) production cross sec- 
tions as well as the top quark branching ratios. The latter are needed for limits from light charged 
Higgs searches in top quark decays. Various input approximations are available, e.g. the production 
cross sections and branching ratios to SM particles can be derived from effective Higgs couplings 
given by the user. HiggsBounds can be called from the command-line, via Fortran90 subrou- 
tines or from a web-interface. It supports^ the SUSY Les Houches Accord (SLHA) for the Minimal 
Supersymmetric Standard Model (MSSM) and the next-to-MSSM (NMSSM). The narrow-width 
approximation is assumed to be applicable. The HiggsBounds result contains information about 
whether the model is excluded at 95% C.L. and which analysis was applied to which Higgs boson. 

The current version is HiggsBounds-3 .8.0 and was released in May 2012. The included 
LHC results are from the 7 TeV run only. The code can be obtained from 

http://higgsbounds. hep forge. org 

A new version HiggsBounds-4 .0.0 including the latest 8 TeV LHC results will be re- 
leased soon. More details will be given in Section 3. 

2. Performance of the Standard Model likeness test 

If a Higgs search is carried out under specific assumptions {e.g. on the CP-properties of the 
Higgs boson or top quark branching ratios) HiggsBounds considers these analyses only if the 
investigated model fulfills these assumptions. In particular, as mentioned in Section 1 , if the Higgs 
seaixh combines several signal topologies under the assumption of the SM, HiggsBounds tests 
whether the model is sufficiently SM-like in these signal topologies. This SM-likeness test is nec- 
essary since the information on efficiencies of the various signal topologies is rarely made publicly 
available. 

For each signal topology /, comprised of the production mode P{h) and decay to the final state 
F, we define an individual signal strength modifier and SM weight w,-, 

[amodel(P(/^))BRn,odel(/i^F)],- ^ ^ [asM{P{H))BRsM{H ^ F)]i 

[asM{P{H))BRsM{H ^ F)]i ' ' Iy[asM(P(//))BRsM(// ^ f )]/ 

Here, the sum runs over all considered signal topologies j of the analysis. Thus, the SM weight ft), 
describes the relative contribution of the signal topology / to the total signal rate in the SM. 

Neglecting the effects from different efficiencies among the signal topologies, the total signal 
strength modifier n is approximated by jJ. = ^, ft),c;. However, this approximation is only valid if 
the signal topologies in the model contribute in similar proportions to the total signal rate as in the 
SM, or in other words, as long as the individual signal strength modifiers c, are approximately the 

^The required HiggsBounds effective coupling input has to be provided via two extra SLHA blocks, see Ref. [4]. 
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Figure 2: Performance of the SM likeness test: We consider the ATLAS H yy search [1] and modify the 
SM normalized squared effective Higgs couplings gli^,^ (left) and gjjyy iV = W,Z) {right) for a Higgs boson 
with mass m = 125 GeV. We show the dependence of the total signal strength modifier jj. and the relevant 
individual signal strength modifiers. The gray regions indicate the parameter space where the SM likeness 
test fails. 



same as the total signal strength modifier /i. Therefore, in order to pass the SM likeness test, we 
require the model to fulfill the following criterion: 



max ft). 



<£, 



with 



and 



:2%, 



(2.2) 



i.e. the maximal weighted deviation of the individual signal strength modifier from the total signal 
strength modifier is required to be less than 2%. We consider e = 2% as a conservative choice, 
considering that the uncertainties on the rate predictions are generally larger. 

An essential part in the calculation of A is the inclusion of the SM weight ft),-, which we want 
to illustrate with an example. We consider again the ATLAS // — )• 77 search [1] and test a model 
with a Higgs boson with a mass m = 125 GeV. We depart from the SM by modifying either 
the (SM normalized) squared effective Higgs coupling to gluons, gjjgg, or to vector bosons, gjivv 
(V = W,Z). All other effective Higgs couplings, in particular the Hjy coupling, are set to their SM 
values. At m = 125 GeV, the SM weights for the LHC at = 7 TeV are 



CO ^ (87.7%, 6.8%, 3.2%, 1.8%, 0.5%) for (ggf, VBF, HW, HZ, Htt) 



(2.3) 



In Fig. 2 we show the signal strength modifier pt and the c, for the signal topologies influenced 
by the modified effective Higgs couplings. Varying gj^^^ influences only the gluon-gluon fusion 
cross section, however, due to its large SM weight, ftjggf « 87.7%, the total signal strength modifier 
pt follows closely c(ggf). The SM likeness test failure at g^^^ = 0.835 and 1.225 is eventually 
caused by the ggf signal topology, although the deviation 5c; for the remaining signal topologies 
VBF, HW , HZ and Htt is much larger here. However, the SM weights of these channels are 
much smaller, thus allowing for a larger deviation. The same effects can be seen when varying 
Shvv ~ W^i-Z)- Now, the c,- of the VBF, HW , HZ signal topologies are affected by the modified 
effective coupling, however, the total signal strength modifier /i is only slightly influenced due the 
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small weight of these channels. Again, the deviation between /i and c(ggf) is eventually causing 
the SM likeness test to fail. Thus, due to the inclusion of the SM weights in Eq. (2.2) subdominant 
signal topologies are allowed to deviate further from jj,. 

The SM weights were introduced with the latest version HiggsBounds-3 .8.0. This led to 
a wider applicability of SM Higgs seaixh results to arbitrary Higgs sectors and thus to a significant 
improvement of the performance of HiggsBounds. 

3. Prospects and new developments 

In July 2012, the long-lasting effort of the Higgs searchers was rewarded with a discovery 
in the Standard Model (SM) Higgs searches for H ^ yy and H — ;> ZZ'*' — 4£ at a mass value 
m « 126 GeV in both LHC experiments ATLAS [5] and CMS [6]. The major tasks after the 
discovery are the determination of the couplings, CP and spin properties of the new Higgs-like 
state. Any deviations from the SM expectation will hint towards new physics beyond the Standai^d 
Model and must be investigated thoroughly. 

The HiggsBounds team is currently developing the program HiggsSignals, aimed to 
perform a test of Higgs sector predictions against the signal rate and mass measurements from 
Higgs collider searches. This information can for instance be used in global fits of models beyond 
the SM or in a generic Higgs coupling determination. In the calculation the systematic un- 
certainties of production cross sections, decay rates, luminosities and Higgs mass predictions are 
taken into account with the full information on their correlations. Furthermore, the code automat- 
ically considers superpositions of the signal rates of Higgs bosons, if their mass difference cannot 
be resolved by the experimental analysis. 

To illustrate some features of the upcoming programs HiggsBounds-4 .0.0 and Higgs- 
Signals-l.O.Owe discuss the interesting possibility of interpreting the LHC Higgs signal as 
the heavier CP-even Higgs boson H of the MSSM (see also Ref. [7, 8, 9] for a detailed analysis). 
We consider the (niA, tanjS) plane for the fixed choice of MSSM parameters 

MsusY = 1 TeV, \Xt\ = 2.4 TeV, n = 1 TeV, Mj = 100 GeV, M2 = 200 GeV, M3 = 800 GeV. 

A definition of these parameters can e.g. be found in Ref. [7]. In Fig. 3 we show the total Ax^ 
likelihood map for this scenario. It is composed of the x^ calculated by HiggsSignals using 
the signal rate and mass measurements presented at the ICHEP2012 conference by ATLAS [5], 
CMS [6], D0 and CDF [10], and a x^ value from the LEP exclusion. This LEP x^ information 
will be included in HiggsBounds-4 .0.0. Furthermore we show the 95% C.L. excluded regions 
obtained with HiggsBounds-3 .8.0. While the ATLAS SM combined analysis of the channels 
H — )■ 77, WW, ZZ [11] and the CMS h,H,A tt search [12] constrain^ the parameter regions 
with tanj8 > 8, the ATLAS charged Higgs search [13] excludes regions at low niA and tanj3. 

We find the best-fit point at {ma, tanj3) = (101.0 GeV, 6.0) with a good fit quality ;t^/ndf = 
29.3/32. It exhibits a charged Higgs boson and a light CP-even Higgs boson h with masses 
mH± = 126 GeV and m/, = 92.3 GeV, respectively. The light Higgs boson h has a reduced coupling 

■'Note that the newer 95% C.L. exclusion limits from the LHC 8 TeV run are even more constraining than the 
exclusion limits applied here. 
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Figure 3: distribution for the heavy Cf -even Higgs boson interpretation of the LHC Higgs signal in the 
(niA, tanjS) plane. The gray striped regions are excluded at 95% C.L. by HiggsBounds. The best-fit point 
is indicated by a green asterisk. 

to Z bosons, hence it escapes the LEP exclusion limits. The rate for the charged Higgs process 
in top decays, t ^ H^b T+V^ft, is fairly close the ATLAS exclusion limit [13]. Thus, with 
more data charged Higgs boson searches will be able to probe this interpretation of the LHC Higgs 
signal. 
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